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Generation of large-scale winds in horizontally anisotropic convection 
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We simulate three-dimensional, horizontally periodic Rayleigh-Benard convection, confined between 
free-slip horizontal plates and rotating about a distant horizontal axis. When both the temperature 
difference between the plates and the rotation rate are sufficiently large, a strong horizontal wind 
is generated that is perpendicular to both the rotation vector and the gravity vector. The wind is 
turbulent, large-scale, and vertically sheared. Horizontal anisotropy, engendered here by rotation, 
appears necessary for such wind generation. Most of the kinetic energy of the flow resides in the 
wind, and the vertical turbulent heat flux is much lower on average than when there is no wind. 


Buoyancy-driven convection is a primary mechanism 
of heat transport in the atmospheres of planets and cool 
stars, the cores of massive stars, and the Earth’s oceans, 
mantle, and outer core. The essential physics of the pro¬ 
cess have most commonly been studied in a layer of fluid 
heated below and cooled above, a configuration called 
Rayleigh-Benard convection (RBC) by Avsec pQ. The 
many natural occurrences of convection and the inherent 
interest of RBC as a complex system have motivated a 
great many laboratory experiments and numerical simu¬ 
lations throughout the past century HE]. 

Large-scale horizontal winds in RBC have often been 
observed, as in the laboratory studies of Malkus [4] and 
Krishnamurti and Howard [5]. The strongest winds yet 
reported occur in simulations of horizontally periodic, 
two-dimensional (2D) RBC between free-slip top and 
bottom boundaries. These winds have a horizontal wave- 
number of zero, and they are vertically sheared, meaning 
that their velocities change with height. For weak ther¬ 
mal forcing, the bifurcations by which ordinary 2D con¬ 
vective cells can lose stability to states with wind have 
been studied using the full governing equations B 13 
and reduced models [SHU]- When the thermal forcing 
is strengthened, the wind in 2D simulations intensifies 
to dominate the underlying convection, and heat trans¬ 
port decreases dramatically [15JI16]. Comparably strong 
winds in 3D RBC have not been reported as yet. 

The simulations described here reveal that strong, 
wavenumber-zero winds can arise also in 3D RBC if the 
horizontal isotropy is broken by uniform rotation about a 
horizontal axis. The resulting flows, which we call windy 
convection, are more complicated than their 2D coun¬ 
terparts but can develop winds of comparable strength. 
Three-dimensional windy convection resembles flow in 
the equatorial regions of rotating spherical shells ([17f 
EDI) and in the outer edge of tokamak plasmas in the 
high-confinement mode m ■ Zonal winds in spherical 
shells are further driven by the way in which rotation 
and curvature conspire to stretch vortices [22], but here 
we find strong winds also without this mechanism. 

Horizontal anisotropy is apparently necessary for 
strong wavenumber-zero winds to persist. The winds 


ceased in our simulations whenever rotation was slowed 
sufficiently, and winds in the initial conditions invariably 
died out in non-rotating simulations [23 . Likewise, so¬ 
lutions with mean winds have been found in truncated 
modal expansions of 3D RBC only when the imposed 
convective pattern is sufficiently anisotropic [24]. Un¬ 
der horizontally isotropic conditions, other sorts of large- 
scale circulation have been seen, but none has had much 
effect on the mean heat flux. Circulation with a hori¬ 
zontal scale comparable to the layer depth can spatially 
organize smaller-scale plumes but does not inhibit their 
heat transport [25] [26], and the wind seen by Krishna¬ 
murti and Howard [5 was apparently not strong enough 
to prevent thermal plumes from traversing the layer. 

We adopt the Boussinesq approximation mm, in 
which the fluid has constant kinematic viscosity, v , ther¬ 
mal diffusivity, x, and coefficient of thermal expansion, 
a. The natural units we use are the layer thickness, d, 
the imposed vertical temperature difference between the 
bounding plates, A, and the thermal time, r t h = d 2 / x. 
Other time scales include the viscous time, r v i s = d 2 /v, 
the dynamical time, Td yn — \/d/ ( gaA ), where g is 
gravitational acceleration, an effective dissipative time, 
Tdiss = y/TthTvis, and the rotation period. 

The ability of buoyancy forces to overcome vis¬ 
cous drag is measured by the Rayleigh number, R = 
i^diss/Tdyn) 2 , and the relative importance of the two dis¬ 
sipative mechanisms is measured by the Prandtl number, 
a = Tfh / T v i s . We set a = 0.71 (air).The dimensionless 
vertical extent is 0 < z < 1, and we set horizontal periods 
of 27r in the x and y directions. The system rotates about 
a distant axis parallel to the y -axis with dimensionless an¬ 
gular speed D, expressed in units of 1 /r t h- The effect of 
rotation is approximated as spatially uniform. 

The dimensionless Boussinesq equations governing the 
velocity u = (' u , v, w), temperature T, and pressure p are 

V • u = 0 (1) 

d t u + u • Vu = — Vp + crV 2 u + crRTz — 2Dy x u (2) 
d t T + u • VT = V 2 T. (3) 

At the top and bottom boundaries, the temperatures are 
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FIG. 1. Instantaneous temperature fields in a vertical slice in the (x, z) plane at (a) t — 0.04 and (b) t — 0.25. The parameters 
are R — 10 7 and 2D = 10 4 , the top and bottom are free-slip, and the rotation vector points into the page. Slices at other values 
of y look qualitatively similar since the rotation makes variation in y weak. 


fixed and uniform, with free-slip velocity conditions: 

T \ z= o = 1/2, T\ z=1 = -1/2 (4) 

du/dz\ z= o,i, dv/dz\ z=0 ,i, w| 2=0 ,i = 0 (5) 

Free-slip boundaries favor strong winds, letting them de¬ 
velop at R that are small enough to be computationally 
accessible. In 2D, winds arise readily even when only one 
boundary is free-slip [29] . 

The governing equations were integrated with a spec¬ 
tral method in the periodic directions with 2/3 dealias¬ 
ing, second-order finite differences in the vertical, and a 
third-order fractional step method in time [23j [30]. We 
used 256 modes in each horizontal direction and 192 ver¬ 
tical levels, spaced more closely near the boundaries. 

We focus first on a simulation with R = 10 7 and 
2D = 10 4 in which a cyclonic wind develops. Cyclonic 
winds, with net vorticity parallel to the rotation vector, 
appear clockwise when observed from the direction of 
negative y. (Anticyclonic winds are also possible and 
are discussed below, along with other rotation rates.) 
Figure 1 shows temperature slices in the (x, z) plane at 
two instants of the simulation. A movie of this simu¬ 
lation is available at [ 51 ] . Beginning from statistically 
homogeneous, random initial conditions, a pair of turbu¬ 
lent convective cells develop, and the temperature field 
is dominated by a hot plume and a cold plume (panel 
a). This stage persists for many turnover times, during 
which the cyclonic cell slowly (and non-monotonically) 
widens while the anticyclonic cell narrows. Eventually, 
the narrowing cell shrinks out of existence as the hot and 
cold plumes collide, and the widening cell breaks through 
to form a cyclonic wind that then dominates the system 
(panel b). Once the wind forms, fluid flows quickly in 
the +x-direction along the top boundary and in the — x- 
direction along the bottom one. Thermal plumes come 


into and out of existence, but they are sheared out hori¬ 
zontally and rarely reach the opposite boundaries. 

The onset of the wind significantly alters not only the 
appearance of the flow but also its integral quantities. 
The specific kinetic energy of motion in the wind direc¬ 
tion, \(u 2 ), is much larger than that in the other two 
directions, |( v 2 + w 2 ). (Angular brackets denote vol¬ 
ume averages.) The decrease in vertical heat flux that 
accompanies the wind’s onset is captured by the Nusselt 
number, N(t) = 1 + (wT), which is the factor by which 
convection amplifies heat transport. 

Figure 2a shows the time variation of the specific ki¬ 
netic energy in different velocity components. In an ini¬ 
tial transient stage the mean flow slowly strengthens. 
The plumes collide near t = 0.08, after which the fully 
developed wind strongly suppresses motions normal to it. 
Still, these motions never cease completely, as shown by 
the small but finite values of v 2 + w 2 ) in Fig. 2a. The 
growth and saturation of the wind are illustrated in Fig. 
2b by vertical profiles of u , averaged over the horizon¬ 
tal directions and 0.025 thermal times. Even during the 
transient stage, the wind becomes fairly strong as the cy¬ 
clonic convective cell grows to dominate the anticyclonic 
one (see Fig. la). Ultimately, the vertical profile of u be¬ 
comes almost linear, with transitory deviations created 
by plume outbursts. 

Figure 2c shows the time evolution of the Nusselt num¬ 
ber for the main experiment (blue line), which first dis¬ 
plays a strong peak as convective plumes emerge explo¬ 
sively from the linear instability. Then N(t) changes lit¬ 
tle while the wind builds in strength. When the wind 
takes over as the cyclonic cell annihilates the anticyclonic 
one, convective heat transport becomes very feeble, as 
reflected by N(t) dropping almost to unity. Heat trans¬ 
port recovers only partially thereafter since the system 
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FIG. 2. (a) Time series of volume-averaged specific kinetic energy in the direction of the wind, |('u 2 }, and normal to the wind, 

|(u 2 + w 2 ); (b) vertical profiles of the wind velocity averaged horizontally and over 0.025 thermal times; (c) time series of the 
Nusselt numbers. An anticyclonic simulation (A) is started at time t = 0.25 by applying the reflection x i —> —x to the cyclonic 
simulation (C). The inset shows mean wind profiles, averaged horizontally and over 0.3 < t < 0.75; (d) vertical profiles of the 
temperature averaged horizontally and over 0.025 thermal times. The dots in panels (a) and (c) mark the starts of the time 
windows over which the profiles in panels (b) and (d) have been averaged. 


remains dominated by the horizontal wind. The time- 
averaged Nusselt number is 7.1 after the wind has fully 
formed (0.20 < t < 0.75), as compared with 22.6 dur¬ 
ing the transient stage (0.05 < t < 0.08). Figure 2d 
shows the evolution of mean vertical temperature pro¬ 
files, which are roughly isothermal in the interior before 
the wind develops but have a strongly unstable stratifi¬ 
cation afterward. 

In addition to being smaller when the wind dominates, 
N(t) varies more strongly; at irregular intervals, there 
are deep minima of convective transport during which 
the Nusselt number dips nearly to its conductive value 
of unity. These minima occur when the wind becomes 
especially strong and nearly shuts off the convection from 
which it draws its energy. Lacking an energy source, the 
wind decays until it is too weak to inhibit the small- 
scale convection, which then rebounds. Such behavior 
is an instance of on-off intermittency [32]. Very strong 
intermittency occurs at similar parameters in 2D RBC 
mm, which may be thought of as an extreme limit of 
horizontal anisotropy in 3D RBC. 

In addition to the cyclonic wind shown in Fig. lb and 
discussed above, there exists at the same parameters a 


second state with sustained anti cyclonic wind — that is, 
with u < 0 near the top and u > 0 near the bottom. We 
found this anticyclonic state by restarting the cyclonic 
simulation at t = 0.25, after reflecting x ha —x. Wind 
profiles and N(t) series for the anticyclonic simulation ac¬ 
company those for the cyclonic one in Fig. 2c. It seems 
there is bistability between the cyclonic and anticyclonic 
attractors, each having its own basin of attraction. Fig. 
2c illustrates that both types of wind are of compara¬ 
ble energy, but the anticyclonic wind bursts much more 
strongly, resembling bursting seen in 2D [15, 16, 291. 

It is rotation that distinguishes between the two wind 
directions. Without rotation, winds in 2D RBC sponta¬ 
neously break the reflectional symmetry x —x to form 
a vertical shear with one sign or the other. Here, Cori¬ 
olis force breaks symmetry not only between the x- and 
//-directions but also between cyclonic and anticyclonic 
winds. For instance, Coriolis forces act upward on fluid 
flowing in the -\-x direction and downward on fluid flow¬ 
ing in the — x direction. (The broken x i-T — x symmetry 
recalls the cyclone-anticyclone asymmetry seen in sta¬ 
bly stratified geophysical turbulence [33].) Coriolis force 
cannot help sustain the wind against viscous dissipation, 
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FIG. 3. (a) vertical profiles of the wind velocity, averaged horizontally and over 0.4 thermal times; (b) ranges of N(t) (mean, 

5% and 95% percentiles, extreme values) explored by the experiments over 0.4 thermal times. Results for rotation rates of 
2Q = 15000, 8000, 5000, 3000,1000, 0 are shown for free-slip boundaries, along with results for 2Q = 10 4 with no-slip (NS) 
boundaries. Cyclonic (dashed lines) and anticyclonic (continuous lines) windy solutions result from different initial conditions. 


however, since it performs no work on the fluid. With 
or without rotation, the wind can be driven only by the 
transfer of energy, via the velocity nonlinearity, from the 
smaller-scale motions that are buoyantly driven. 

Keeping R = 10 7 and free-slip boundaries, we have 
explored the effects of rotation rate. Figure 3a gives 
mean horizontal velocity profiles for both cyclonic and 
anticyclonic states at various rotation rates. Bistability 
between each type of wind was found in all cases with 
2Q > 3000. (The windy state of Fig. lb was used to 
initialize the cyclonic simulations, and its x — x re¬ 
flection was used to initialize most anticyclonic ones.) In 
both the cyclonic and anticyclonic states, the mean winds 
grow as rotation is strengthened. Fully developed windy 
convection like that in Fig. lb, where the wind shear is 
strong enough to inhibit convection rolls, occurs only for 
sufficiently high rotation rates. In our simulations, true 
windy convection occurs only when 2Q > 5000 in the 
anticyclonic cases and when 2D > 8000 in the cyclonic 
ones. At lower rotation rates we always find a pair of 
convective rolls, coexisting with a mean wind that is too 
weak to destroy them. Fig. 3a also reports a simulation 
with no-slip boundary conditions. A slight wind develops 
in this case, but it is not nearly strong enough to destroy 
the convective cells and create true windy convection. 

Figure 3b summarizes the ranges of N(t) found in sim¬ 
ulations at various 2D. At the lower rotation rates where 
rolls persist, N(t) behaves much as it does in the non¬ 
rotating case, where there is no mean wind. When the 
rotation rate is high enough for the wind to preclude con¬ 
vective overturning, we find the dramatically altered heat 
transport described above for our main simulation: N(t) 
is highly intermittent, and its time average is small. It 
seems that the wind is strong enough to significantly alter 
N(t) only when it is strong enough to prohibit rolls—that 
is, when u has a single sign near the top boundary and 
the opposite sign near the bottom one. 


In summary, we have observed the formation of strong 
horizontal winds in 3D RBC when a forced symmetry 
breaking, here brought about by imposed rotation, is 
strong enough. When rotation is too weak for wind to 
be sustained, any wind that starts to form is quickly dis¬ 
rupted by motions normal to it. This is natural because 
the wind does not suppress convective instability normal 
to it, whereas the wind does partially suppress motions 
that are not normal to it and from which it draws its 
energy. This effect has been demonstrated in the related 
system where Couette shear is imposed on RBC mm- 
The tendency of strong winds to quench the motions that 
generate them has analogs also in the large-scale mag¬ 
netic fields of convective dynamo theory [36, 37, and the 
zonal flows of tokamak plasmas [38] . 

Our findings support the conjecture that the upscale 
cascade of energy in anisotropic turbulent convection, 
which here drives sheared winds, drives differential ro¬ 
tation in the equatorial regions of planetary atmospheres 
and stellar convective zones. Such differential rotation 
may have interesting consequences. At the top and bot¬ 
tom of the solar convection zone, for example, shear flows 
have been detected. The shearing flow at the bottom - 
the tachocline [39 - is suspected of generating the solar 
activity that is manifested by sunspots [40 , and the in¬ 
termittent behavior of the heat transport described here 
recalls the irregularity of such activity m- Other pos¬ 
sible locations of windy convection are in the convective 
cores of massive stars that rotate rapidly. Large-scale 
shear arising there might influence nuclear reaction rates 
[42] . with implications for stellar evolution. More im¬ 
mediately, the desire to extrapolate windy convection to 
astrophysically or geophysically interesting regimes calls 
for further exploration of parameter space. 
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